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Abstract 

We recall that, at both the intermittency transitions and at the 
Feigenbaum attractor in unimodal maps of non-linearity of order £ > 
1, the dynamics rigorously obeys the Tsallis statistics. We account for 
the q-indices and the generalized Lyapunov coefficients X q that charac- 
terize the universality classes of the pitchfork and tangent bifurcations. 
We identify the Mori singularities in the Lyapunov spectrum at the 
edge of chaos with the appearance of a special value for the entropic 
index q. The physical area of the Tsallis statistics is further probed 
by considering the dynamics near criticality and glass formation in 
thermal systems. In both cases a close connection is made with states 
in unimodal maps with vanishing Lyapunov coefficients. 
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1 Introduction 

Ever since the proposition, in 1988, by Tsallis [1], [2] of a nonextensive gener- 
alization of the canonical formalism of statistical mechanics a spirited discus- 
sion [3], [4] has grown on the foundations of this branch of physics. But, only 
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until recently [5]- [12] has there appeared firm evidence about the relevance of 
the generalized formalism for specific model system situations. These studies 
present rigorous results on the dynamics associated to critical attractors in 
prototypical nonlinear one-dimensional maps, such as those at the pitchfork 
and tangent bifurcations and at the accumulation point of the former, the 
so-called onset of chaos [13]. As these are very familiar and well understood 
features of these maps it is of interest to see how previous knowledge fits 
in with the new perspective. Also, clear-cut calculations may help clarify 
the physical reasons (believed in these examples to be a breakdown in the 
chain of increasing randomness from non-ergodicity to completely developed 
chaoticity) for the departure from the Boltzmann-Gibbs (BG) statistics and 
the competence of the non-extensive generalization. Here we review briefly 
specific results associated to the aforementioned states in unimodal maps 
which are all characterized by vanishing Lyapunov coefficients. We also re- 
call how the dynamics at the tangent bifurcation appears to be related to that 
at thermal critical states. In addition, time evolution at the noise-perturbed 
onset of chaos is seen to be closely analogous to the glassy dynamics observed 
in supercooled molecular liquids. 

2 Dynamics at the pitchfork and tangent bi- 
furcations 

As pointed out in Refs. [5], [6], [9] the long-known [13] exact geometric or 
static solution of the Feigenbaum renormalization group (RG) equations for 
the tangent bifurcation in unimodal maps of nonlinearity £ > 1 also de- 
scribes the dynamics of iterates at such state. A straightforward extension 
of this approach applies also to the pitchfork bifurcations. We recall that 
the period-doubling and intermittency transitions are based on the pitchfork 
and the tangent bifurcations, respectively, and that at these critical states the 
ordinary Lyapunov coefficient Ai vanishes. The sensitivity to initial condi- 
tions £ t was determined analytically and its relation with the rate of entropy 
production examined [9]. The fixed-point expressions have the specific form 
that corresponds to the temporal evolution suggested by the nonextensive 
formalism. These studies contain the derivation of their g-generalized Lya- 
punov coefficients \ q and the description of the different possible types of 



2 



sensitivity £ t . 

By considering as starting point the (- logistic map f^x) = 1 — /i [xf, ( > 
1, — 1 < x < 1, it is found that for both infinite sets of pitchfork and tangent 
bifurcations £ t , defined as £ t = lim/± X0 ^ (Ax t / Ax ) (where Ax is the length 
of the initial interval and Ax t its length at time t), has the form suggested 
by Tsallis, 

&(xo) = exp ? [A ? (x ) t] = [l-(q- l)X q (x ) t]~^, (1) 

that yields the customary exponential £ t with Lyapunov coefficient Xi(x ) 
when 5 — > 1. In Eq. (1) q is the entropic index and X q is the g-generalized 
Lyapunov coefficient; expire) = [1 — (q — l)x]~ 1 ^ q ~ 1 ' ) is the g-exponential 
function. The pitchfork and the left-hand side of the tangent bifurcations 
display weak insensitivity to initial condition, while the right-hand side of 
the tangent bifurcations presents a 'super-strong' (faster than exponential) 
sensitivity to initial conditions [6]. 

For the transition to periodicity of order n the composition ffp is first 
considered. In the neighborhood of one of the n points tangent to the line 
with unit slope one obtains f( n \x) = x + u \x\ z + o(|x| 2 ), where u > is the 
expansion coefficient. The general result obtained is q = 2 — and X q (xo) = 
zuXq~ 1 [Q\, [9]. At the tangent bifurcations one has f^(x) = x + ux 2 + 
o(x 2 ), u > 0, and from z = 2 one gets q = 3/2. For the pitchfork bifurcations 
one has instead fffi(x) = x + ux 3 + o(x 3 ), because cPff 2 ) /dx 2 = at these 
transitions, and u < is now the coefficient associated to d 3 fj 2 ^ / dx 3 < 0. 
In this case we have z = 3inq = 2 — z^ 1 and one obtains q = 5/3. Notably, 
these specific results for the index q are valid for all £ > 1 and therefore 
define the existence of only two universality classes for unimodal maps, one 
for the tangent and the other one for the pitchfork bifurcations [6]. See Figs. 
2 and 3 in Ref. [6]. 

Notice that our treatment of the tangent bifurcation differs from other 
studies of intermittency transitions [14] in that there is no feed back mecha- 
nism of iterates into the origin of f^(x) or of its associated fixed-point map. 
Here impeded or incomplete mixing in phase space (a small interval neigh- 
borhood around x — 0) arises from the special 'tangency' shape of the map at 
the pitchfork and tangent transitions that produces monotonic trajectories. 
This has the effect of confining or expelling trajectories causing anomalous 
phase-space sampling, in contrast to the thorough coverage in generic states 
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with Ai > 0. By construction the dynamics at the intermittency transitions, 
describe a purely nonextensive regime. 



3 Link between intermittent dynamics and 
dynamics at thermal criticality 

An unanticipated relationship has been shown by Contoyiannis et al [15]- [18] 
to exist between the intermittent dynamics of nonlinear maps at a tangent 
bifurcation and the equilibrium dynamics of fluctuations at an ordinary ther- 
mal critical point. With the aim of obtaining the properties of clusters or 
domains of the order parameter at criticality, a Landau- Ginzburg- Wilson 
(LGW) approach was carried out such that the dominant contributions to 
the partition function arise from a singularity (similar to an instanton) lo- 
cated in the space outside the cluster [15], [16]. Then [17], [18], a nonlinear 
map for the average order parameter was constructed whose dynamics repro- 
duce the averages of the cluster critical properties. This map has as a main 
feature the tangent bifurcation and as a result time evolution is intermittent. 

The starting point is the partition function of the ci-dimensional system 
at criticality, 

Z = J L>[0]exp(-r c [0]), (2) 

where 

r c [0]=si J Q dv[l(V<t>) 2 + g 2 \<l>\ s+1 ] (3) 

is the critical LGW free energy of a system of (/-dimensional volume Q, <fi 
is the order parameter (e.g. magnetization per unit volume) and 5 is the 
critical isotherm exponent. By considering the space-averaged magnetization 
$ = J v <f)(x)dV, the statistical weight 

p($) = exp(-r c [S])/Z, (4) 

where r c [$] ~ g\g2^ s+1 and Z = J c?$ exp(— r c [$]), was shown to be the in- 
variant density of a statistically equivalent one-dimensional map. The func- 
tional form of this map was obtained as the solution of an inverse Frobenius- 
Perron problem [17]. For small values of $ the map has the form 

<fn+l = <fn + M< +1 +e, (5) 
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where the amplitude u depends on gi, (^and 5, and the shift parameter 
e ~ R~ d . Eq. (5) can be recognized as that describing the intermittency 
route to chaos in the vicinity of a tangent bifurcation [13]. The complete 
form of the map displays a superexponentially decreasing region that takes 
back the iterate close to the origin in one step. Thus the parameters of 
the thermal system determine the dynamics of the map. Averages made 
of order-parameter critical configurations are equivalent to iteration time 
averages along the trajectories of the map close to the tangent bifurcation. 
The mean number of iterations in the laminar region was seen to be related 
to the mean magnetization within a critical cluster of radius R. There is a 
corresponding power law dependence of the duration of the laminar region 
on the shift parameter e of the map [17]. For e > the (small) Lyapunov 
coefficient is simply related to the critical exponent 5 [16]. 

As the size of subsystems or domains of the critical system is allowed to 
become infinitely large the Lyapunov coefficient vanishes and the duration of 
the laminar episodes of intermittency diverge. Without the feedback to the 
origin feature in the map we recover the conditions described in the previous 
section. See Ref. [10] for more details. 

4 Dynamics inside the Feigenbaum attractor 

The dynamics at the chaos threshold, also referred to as the Feigenbaum 
attractor, of the ^-logistic map at fi c has been analyzed recently [7], [8]. By 
taking as initial condition xo = we found that the resulting orbit consists of 
trajectories made of intertwined power laws that asymptotically reproduce 
the entire period-doubling cascade that occurs for \i < ji c . This orbit captures 
the properties of the so-called 'superstable' orbits at ~p, n < /i c , n — 1,2,... 
[13]. Here again the Lyapunov coefficient Ai vanishes and in its place there 
appears a spectrum of g-Lyapunov coefficients X^ k \ This spectrum was orig- 
inally studied in Refs. [20], [21] and our interest has been to examine the 
relationship of its properties with the Tsallis statistics. We found that the 
sensitivity to initial conditions has precisely the form of a g-exponential, of 
which we determine the g-index and the associated A^. The appearance of 
a specific value for the q index (and actually also that for its conjugate value 
Q = 2 — q) turns out to be due to the occurrence of Mori's 'g transitions' [21] 
between 'local attractor structures' at ji c . Furthermore, we have also shown 
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that the dynamical and entropic properties at // c are naturally linked through 
the nonextensive expressions for the sensitivity to initial conditions £j and 
for the entropy S q in the rate of entropy production K q k \ We have corrob- 
orated analytically [7], the equality A^ = given by the nonextensive 
statistical mechanics. Our results support the validity of the g-generalized 
Pesin identity at critical points of unimodal maps. 

Thus, the absolute values for the positions x T of the trajectory with x t =o = 
at time-shifted r = t + 1 have a structure consisting of subsequences with a 
common power- law decay of the formr -1 / 1-9 with q — 1 — In 2/(£ — 1) In a(C), 
where is the Feigenbaum universal constant that measures the period- 
doubling amplification of iterate positions [6]. That is, the Feigenbaum 
attractor can be decomposed into position subsequences generated by the 
time subsequences r = (2k + 1)2™, each obtained by proceeding through 
n — 0, 1, 2, ... for a fixed value of k = 0, 1, 2, .... See Fig. 1. The k — subse- 
quence can be written as %t = exp 2 _ ? (— X q °H) with \ q ^ = (£ — 1) \na(Q/ In 2. 
These properties follow from the use of xo = in the scaling relation [6] 

(6) 

The sensitivity associated to trajectories with other starting points x ^ 
within the attractor can be determined similarly with the use of the time 
subsequences r = (2A;+l)2 n . One obtains \ q k) = ((-1) lna(C)/(2fc+l) In 2 > 
0, k — 0,1,2,..., the positive branch of the Lyapunov spectrum, when the 
trajectories start at the most crowded (x t= q = 1) and finish at the most 
sparse (x T= 2" = 0) region of the attractor. By inverting the situation we 
obtain \^} q = -2(£ - l)lna/(2A; + l)ln2 < 0, k = 0,1,2,..., the negative 
branch of \ q k \ i.e. starting at the most sparse (x T= o = 0) and finishing at the 
most crowded (x T=2 n+i = 0) region of the attractor. Notice that Q = 2 — q as 
expg(y) = 1/ exp g (— y). For the case ( = 2 see Refs. [7] and [8], for general 
C > 1 see Refs. [11] and [19] where also a different and more direct derivation 
is used. So, when considering these two dominant families of orbits all the 
g-Lyapunov coefficients appear associated to only two specific values of the 
Tsallis index, q and 2 — q. 

As a function of the running variable — oo < q < oo the coefficients 
become a function A(q) with two steps located at q = g = 1 =F ln2/(^ — 
l)lna(£). In this manner contact can be established with the formalism 
developed by Mori and coworkers and the q phase transition obtained in Ref. 
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Figure 1: Absolute values of positions in logarithmic scales of the first 1000 
iterations r for a trajectory of the logistic map at the onset of chaos /x c (0) 
with initial condition x in = 0. The numbers correspond to iteration times. 
The power-law decay of the time subsequences described in the text can be 
clearly appreciated. 

[21]. The step function for A(q) can be integrated to obtain the spectrum 
0(q) (A(q) = d(j)/d\(q)) and its Legendre transform tp(X) (= — (1 — q)A), 
the dynamic counterparts of the Renyi dimensions D q and the spectrum f(a) 
that characterize the geometry of the attractor. The constant slopes in the 
spectrum represent the Mori's q transitions for this attractor and the 

value 1 — q coincides with that of the slope previously detected [20], [21]. See 
Fig. 2. Details appear in Ref. [19]. 

Ensembles of trajectories with starting points close to x T= o = 1 expand in 
such a way that a uniform distribution of initial conditions remains uniform 
for all later times t < T where T marks the crossover to an asymptotic regime. 
As a consequence of this we established [8] the identity of the rate of entropy 
production with A^. The g-generalized rate of entropy production K q 
is defined via K q t = S q (t) — 5^(0), t large, where 

is the Tsallis entropy, and where ln q y = (y 1 ^ — 1)/(1 — q) is the inverse of 
exp q (y). See Figs. 2 and 3 in Ref. [8]. 
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Figure 2: a) The Lyapunov coefficient function A(q) at the chaos threshold 
at fi c and b) the spectrum ^(A). See text for description 

Consider now the logistic map ( = 2 in the presence of additive noise 

x t +i = U(x t ) = 1 - fJ^t + Xt&, -1 < x t < 1, < n < 2, (8) 

where Xt is Gaussian-distributed with average (xtXf) — $t.t', and a measures 
the noise intensity. We recall briefly the known properties of this problem 
[13], [22]. Except for a set of zero measure, all the trajectories with fi c (a = 
0) and initial condition — 1 < xq < 1 fall into the attractor with fractal 
dimension df — 0.5338.... These trajectories represent nonergodic states, 
since as t — > oo only a Cantor set of positions is accessible out of the total 
phase space — 1 < x < 1. For o > the noise fluctuations wipe the sharp 
features of the periodic attractors as these widen into bands similar to those 
in the chaotic attractors, nevertheless there remains a well-defined transition 
to chaos at // c (<t) where the Lyapunov exponent Ai changes sign. The period 
doubling of bands ends at a finite value 2 N ^ as the edge of chaos transition is 
approached and then decreases at the other side of the transition. This effect 
displays scaling features and is referred to as the bifurcation gap [13], [22]. 
When a > the trajectories visit sequentially a set of 2 n disjoint bands or 
segments leading to a cycle, but the behavior inside each band is completely 
chaotic. These trajectories represent ergodic states as the accessible positions 
have a fractal dimension equal to the dimension of phase space. Thus the 
removal of the noise a — > leads to an ergodic to nonergodic transition in 
the map. 
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In the presence of noise (a small) one obtains instead of Eq. (6) [12] 

g ( T Vi-« x ) + xar^G^r 1 / 1 '^) , (9) 



X T = T- 1 ^ 



or 



where G\(x) is the first order perturbation eigenfunction, and where r 
1 - In 2/ In k ~ 0.6332. Use of x = yields x T = r^ 1 / 1 ^ 1 + x^ 1/1_r 
x t = exp 2 _ q (—\ q t) [1 + xo" exp r (A r t)] where t = r — 1 and A r = In In 2. 
At each noise level a there is a 'crossover' or 'relaxation' time t x — t x — 1 
when the fluctuations start suppressing the fine structure of the orbits with 
Xq = 0. This time is given by t x = cr r ~ 1 , the time when the fluctuation 
term in the perturbation expression for x T becomes unbounded by a, i.e. 
x Tx = r~ 1 / 1 ~ q |1 + x\- There are two regimes for time evolution at n c {a). 
When t < t x the fluctuations are smaller than the distances between neigh- 
boring subsequence positions of the a = orbit at fJ> c (0), and the iterate 
positions with o > fall within small non overlapping bands each around 
the a = position for that r. Time evolution follows a subsequence pat- 
tern analogous to that in the noiseless case. When r ~ r x the width of the 
noise-generated band reached at time t x = 2 N matches the distance between 
adjacent positions where iV ~ — mcr/mre, and this implies a cutoff in the 
progress along the position subsequences. At longer times r > t x the orbits 
no longer follow the detailed period-doubling structure of the attractor. The 
iterates now trail through increasingly chaotic trajectories as bands merge 
with time. This is the dynamical image - observed along the time evolution 
for the orbits of a single state /x c (°") - of the static bifurcation gap originally 
described in terms of the variation of the control parameter \x [22], [23], [24]. 
The plateau structure of relaxation and the crossover time t x can be clearly 
observed in Fig. lb in Ref. [25] where < x\ > — < x t > 2 is shown for several 
values of a. 



5 Parallels with dynamics near glass forma- 
tion 

We recall the main dynamical properties displayed by supercooled liquids on 
approach to glass formation. One is the growth of a plateau and for that 
reason a two-step process of relaxation in the time evolution of two-time cor- 
relations [26]. This consists of a primary power-law decay in time difference 
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At (so-called (3 relaxation) that leads into the plateau, the duration t x of 
which diverges also as a power law of the difference T — T g as the temper- 
ature T decreases to a glass temperature T g . After t x there is a secondary 
power law decay (so-called a relaxation) away from the plateau [26]. A sec- 
ond important (nonequilibrium) dynamic property of glasses is the loss of 
time translation invariance observed for T below T g , a characteristic known 
as aging [27]. The time fall off of relaxation functions and correlations display 
a scaling dependence on the ratio t/t w where t w is a waiting time. A third 
notable property is that the experimentally observed relaxation behavior of 
supercooled liquids is effectively described, via reasonable heat capacity as- 
sumptions [26], by the so-called Adam-Gibbs equation, t x = Aexp(B /TS C ), 
where t x is the relaxation time at T, and the configurational entropy S c is 
related to the number of minima of the fluid's potential energy surface [26]. 
We compare the dynamic properties at the edge of chaos described in the 
previous section with those known for the process of vitrification of a liquid 
as T — > T g . 

At noise level a the orbits visit points within the set of 2 N bands and, as 
explained in Ref. [12], this takes place in time in the same way that period 
doubling and band merging proceeds in the presence of a bifurcation gap 
when the control parameter is run through the interval < fi < 2. Namely, 
the trajectories starting at xq = duplicate the number of visited bands at 
times t = 2 n , n = 1, N, the bifurcation gap is reached at r x = 2 N , after 
which the orbits fall within bands that merge by pairs at times r = 2 N+n , 
n = 1,...,N. The sensitivity to initial conditions grows as £ t = exp q (\ q t) 
(q = 1 — In 2/ In a < 1) for t < t x , but for t > t x the fluctuations dominate 
and £t grows exponentially as the trajectory has become chaotic (q = 1) 
[12]. This behavior was interpreted [12] to be the dynamical system analog 
of the a relaxation in supercooled fluids. The plateau duration t x — > oo as 
a — > 0. Additionally, trajectories with initial conditions x not belonging to 
the attractor exhibit an initial relaxation process towards the plateau as the 
orbit approaches the attractor. This is the map analog of the f3 relaxation 
in supercooled liquids. 

The entropy S c (fi c (a)) associated to the distribution of iterate positions 
(configurations) within the set of 2 N bands was determined in Ref. [12]. 
This entropy has the form S' c (/i c (o r )) = 2 N as, since each of the 2^ bands 
contributes with an entropy as, where s = — j\p(x) m P(x)^X an d where 
p(x) is the distribution for the noise random variable. Given that 2 N = l+t x 
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and a = (1 + t x ) 1 / 1 r , one has S c (/j, c , t x )/s — (1 + t x ) r ^ r or, conversely, 

t x = ( S /^c) (1 - r)/r - (10) 

Since t x ~ o" r ~\ r — 1 ~ —0.3668 and (1 — r)/r ~ 0.5792 then ^ — > oo and 
S c — > as <y — > 0, i.e. the relaxation time diverges as the 'landscape' entropy 
vanishes. We interpret this relationship between t x and the entropy S c to be 
the dynamical system analog of the Adam-Gibbs formula for a supercooled 
liquid. Notice that Eq.(10) is a power law in S" 1 while for structural glasses 
it is an exponential in S*" 1 [26]. This difference is significant as it indicates 
how the superposition of molecular structure and dynamics upon the bare 
ergodicity breakdown phenomenon built in the map modifies the vitrification 
properties. 

The aging scaling property of the trajectories Xt at /i c (cr) was examined 
in Ref. [12]. The case a = is readily understood because this property is 
actually built into the position subsequences x T = <?^(0) , r = (2k + l)2 n , 
k, n — 0, 1, ... referred to above. These subsequences can be employed for the 
description of trajectories that are at first held at a given attractor position 
for a waiting period of time t w and then released to the normal iterative 
procedure. For illustrative purposes we select the holding positions to be 
any of those for a waiting time t w — 2k + 1, k — 0, 1, ... and notice that for 
the x in = orbit these positions are visited at odd iteration times. The lower- 
bound positions for these trajectories are given by those of the subsequences 
at times (2k + 1)2™. See Fig. 1. Writing r as r = t w + t we have that 
t/t w = 2 n -l and x t+tw = g^(0)g^(0) or 

x t+tw =g^\0)exp q (-\ q t/t w ). (11) 

This fully developed aging property is gradually modified when noise is 
turned on. The presence of a bifurcation gap limits its range of validity to 
total times t w + 1 < t x (a) and so progressively disappears as a is increased. 



6 Concluding remarks 

The implications of joining the results about intermittency described in Sec- 
tions 2 and 3 are apparent. In the critical clusters of infinite size R — * oo the 
dynamics of fluctuations obeys the nonextensive statistics. This is expressed 
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via the time series of the average order parameter $ n , i.e. trajectories $ n 
with close initial values separate in a superexponential fashion according to 
Eq. (1) with q = (25 + 1)/(S + 1) > 1 and with a g-Lyapunov coefficient X q 
determined by the system parameter values 5, g±, g<i [10]. 

Also, as described in Sections 4 and 5, the dynamics of noise-perturbed 
logistic maps at the chaos threshold exhibits the most prominent features 
of glassy dynamics in supercooled liquids. The existence of this analogy 
cannot be considered accidental since the limit of vanishing noise amplitude 
(j^O (the counterpart of the limit T — T g — > in the supercooled liquid) 
entails loss of ergodicity. Here we have shown that this nonergodic state 
corresponds to the limiting state, a — > 0, t x — > oo, for a family of small a 
noisy states with glassy properties, that are noticeably described for t < t x 
via the g-exponentials of the nonextensive formalism [12]. 

What is the significance of the connections we have reviewed? Are there 
other connections between critical phenomena and transitions to chaos? Are 
all critical states - infinite correlation length with vanishing Lyapunov coef- 
ficients - outside BG statistics? Where, and in that case why, does Tsallis 
statistics apply? Is ergodicity failure the basic playground for applicability 
of generalized statistics? We can mention some noticeable limitations in the 
examples discussed. In the case of intermittency we have focused only on 
a single laminar episode (or scape from a position very close to tangency), 
though this can be of very large duration. In the case of dynamics at the edge 
of chaos we have payed attention only to the dominant features of the mul- 
tifractal attractor (the most sparse x = and most crowded x = 1 regions) 
as starting and final orbit positions. 
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